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A TUBE FORMULA FOR THE KOCH SNOWFLAKE CURVE, 
WITH APPLICATIONS TO COMPLEX DIMENSIONS. 

MICHEL L. LAPIDUS AND ERIN P. J. PEARSE 



Abstract 

A formula for the interior e-neighbourhood of the classical von Koch snowflake curve is computed 
in detail. This function of e is shown to match quite closely with earlier predictions from [La-vFl] 
of what it should be, but is also much more precise. The resulting 'tube formula' is expressed 
in terms of the Fourier coefficients of a suitable nonlinear and periodic analogue of the standard 
Cantor staircase function and reflects the self-similarity of the Koch curve. As a consequence, the 
possible complex dimensions of the Koch snowflake are computed explicitly. 



1. Introduction 

In [La-vFl], the authors lay the foundations for a theory of complex dimensions 
with a rather thorough investigation of the theory of fraetal strings (see also, e.g., 
[LaPol-2,LaMa,Lal-3,IIaLa,IIeLa,La-vF2-3]); that is, fractal subsets of R. Such 
an object may be represented by a sequence of bounded open intervals of length Ij : 

oo 

£:={lj}T=i: with^Z,<oo. (1.1) 

J=l 

The authors are able to relate geometric and physical properties of such objects 
through the use of zeta functions which contain geometric and spectral information 
about the given string. This information includes the dimension and measurability 
of the fractal under consideration, which we now recall. 

For a nonempty bounded open set 17 C M, V{e) is defined to be the inner 
e -neighborhood of O: 

V{e) ■.= yoli{x en : d{x,dn) < e}, (1.2) 

where voli denotes 1-dimensional Lebesgue measure. Then the Minkowski dimen- 
sion of the boundary dil (i.e., of the fractal string £) is 

D = Dan ^ nif{t > : V{e) = O (e^^*) as e ^ 0+}. (1.3) 

Finally, dQ is Minkowski measurable if and only if the limit 

M = M{D; dn) = lin-Y V(£)e"(^"-°) (1.4) 

exists, and lies in (0, oo). In this case, A4 is called the Minkowski content of Sfi. 



2000 Mathematics Subject Classification 26A30, 28A12, 28A75, 28A80 (primary); 11K55, 26A27, 
28A78, 28D20, 42A16 (secondary). 

The work of MLL was partially supported by the US National Science Foundation under grant 
DMS-0070497. Support by the Centre Emile Borel of the Institut Henri Poincare (IHP) in Paris 
and by the Institut des Hautes Etudes Scientiflques (IHES) in Bures-sur-Yvette, France, during 
completion of this article is also gratefully acknowledged by MLL. 



2 



M. L. LAPIDUS AND E. P. J. PEARSE 



If, more generally, fl is an open subset of M.'^, then analogous definitions hold 
if 1 is replaced by d in (1.2)-(1.4). Hence, void denotes the d-dimcnsional volume 
(which is area for d = 2) in the counterpart of (1.2) and e^~*, e^^^^-'^' are replaced 
by s'^-* in (1.3), and e-('^-^) in (1.4), respectively. In (1.2)-(1.4), d = 1 and the 
positive numbers Ij are the lengths of the connected components (open intervals) of 
fl, written in nonincreasing order. In much of the rest of the paper (where fl is the 
snowflake domain of M^), we have d = 2. See, e.g., [Man, Tr, Lai, LaPol-2,Mat,La- 
vFl] and the relevant references therein for further information on the notions 
of Minkowski-Bouligand dimension (also called 'box dimension') and Minkowski 
content. 

The conipk;x dimensions of a fractal string C are defined to be the poles (of the 
meromorphic continuation) of its geometric zeta function 

oo 

Cc{s) = Y.lh (1.5) 

in accordance with the result that 

oo 

Dc = inf{c7 > : < oo}, (1.6) 

i.e., that the Minkowski dimension of a fractal string is the abscissa of convergence 

of its geometric zeta function [La2] . 

To rephrase, we define the complex dimensions of C to be the set 

I? = {w e C : is defined and has a pole at u)}. (1.7) 

One reason why these complex dimensions are important is the (explicit) tubular 
formula for fractal strings, a key result of [La-vFl]. Namely, that under suitable 
conditions on the string £, we have the following tube formula: 

^(-) = E--|S+^(-)' (1-8) 

where the sum is taken over the complex dimensions w of £, and the error term 
R(e) is of lower order than the sum as e ^ 0^. (Sec [La-vFl, Thm. 6.1, p. 144].) 
In the case when £ is a self-similar string i.e., when dL is a self-similar subset of 
K, one distinguishes two complementary cases (see [L,La3] and [La-vFl, §2.3] for 
further discussion of the lattice/nonlattice dichotomy): 

(1) In the lattice case, i.e., when the imdcrlying scaling ratios are rationally 
dependent, the error term is shown to vanish identically and the complex 
dimensions lie periodically on vertical lines (including the line Res = D). 

(2) In the nonlattice case, the complex dimensions are quasiperiodically dis- 
tributed and s = £> is the only complex dimension with real part D. Es- 
timates for R{£) are given in [La-vFl, Thm. 6.20, p. 154] and more precisely 
in [La-vF2-3]. Also, C is Minkowski measurable if and only if it is nonlat- 
tice. See [La-vFl, Chap. 2 and Chap. 6] for details, including a discussion of 
quasiperiodicity. 

These results pertain only to fractal subsets of M, although since this paper was 
completed and using a different approach, some preliminary work on the extension 
to the higher-dimensional case has been done in [Pe] and [LaPel-2]. The founda- 
tional special case of 'fractal sprays' [LaPo2] is discussed in [La-vFl, §1.4]. 




Figure 1. The Koch curve K and Koch snowflake domain Q.. 



It is the aim of the present paper to make some first steps in this direction. We 
compute V{e) for a well-known (and well-studied) example, the Koch snowflake, 
with the hope that it may help in the development of a general higher-dimensional 
theory of complex dimensions. This curve provides an example of a lattice self- 
similar fractal and a nowhere differentiable plane curve. Further, the Koch snowflake 
can be viewed as the boundary dVl of a bounded and simply connected open set 
J7 C and that it is obtained by fltting together three congruent copies of the 
Koch curve as shown in Fig. 1. A general discussion of the Koch curve may be 
found in [Man, §11.6] or [Fa, Intro, and Chap. 9]. 

The Koch curve is a self-similar fractal with dimension D logg 4 (Hausdorff 
and Minkowski dimensions coincide for the Koch curve) and may be constructed 
by means of its self-similar structure (as in [Ki, p. 15] or [Fa, Chap. 9]) as follows: 
let Q — ^-\- with i — \/— T, and define two maps on C by 

fi{z):=QZ and f2{z) := {I ~ g){z ~ I) + I. 

Then the Koch curve is the self-similar set of with respect to {/i, /2}; i.e., the 
unique nonempty compact set ii' C satisfying K = fi {K) U /2 {K) . 
In this paper, we prove the following new result: 

Theorem 1.1. The area of the inner e -neighbourhood of the Koch snowflake is 
given by the following tube formula: 

T/(£) = Gl(£)£2-^ + G'2(£)e^ (1.9) 

where D — logg 4 is the Minkowski dimension of d^l, p := 27r/ log 3 is the oscillatory 
period, and Gi and G2 are periodic functions (of multiplicative period 3) which are 
discussed in full detail in Thm. 5.1. This formula may also be written 

V{e) = + V-n^'""^, (1-10) 

for suitable constants V'n which depend only on n. These constants are expressed 
in terms of the Fourier coefficients g^ of a multiplicative function which bears struc- 
tural similarities to the classical Cantor-Lebesgue function described in more detail 
in §6. 

While this formula is new, it should be noted that a previous approximation 
has been obtained in [La-vFl, §10.3]; see (2.10) in Rem. 2.3 below. Our present 
formula, however, is exact. By reading off the powers of e appearing in (1.10), we 
immediately obtain the following corollary: 
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Corollary 1.2. The possible complex dimensions of the Koch snowflake are 



This is illustrated in Fig. 8. Also, for more precision regarding Cor. 1.2, see 
Rem. 5.3 below, as well as the discussion surrounding (5.5). 

Remark 1.3. The significance of the tube formula (1.9) is that it gives a de- 
tailed account of the oscillations that are intrinsic to the geometry of the Koch 
snowflake curve. More precisely, the real part D yields the order of the amplitude 
of these oscillations (as a function of e) while the imaginary part np = 27rri/log3 
{n = 0,1,2,...) gives their frequencies. This is in agreement with the 'philoso- 
phy' of the mathematical theory of the complex dimensions of fractal strings as 
developed in [La-vFl]. Additionally, if one can show the existence of a complex 
dimension with real part D and imaginary part inp, n ^ 0, then Theorem 1.1 im- 
mediately implies that the Koch curve is not Minkowski measurable, as conjectured 
in [La3, Conj. 2&3, pp. 159,163-4]. 

The rest of this paper is dedicated to the proof of Thm. 1.1 (stated more precisely 
as Thm. 5.1). More specifically, in §2 we approximate the area V{e) of the inner 
e-neighborhood. In §3 we take into account the 'error' resulting from this approx- 
imation. We study the form of this error in §3.1, and the amount of it in §3.2. In 
§4 we combine §2 and §3 to deduce the tube formula (1.9) in the more precise form 
given in §5. We also include in §5 some comments on the interpretation of Thm. 5.1. 
Finally, in §6 we sketch the graph and briefly discuss some of the properties of the 
Cantor-like and multiplicatively periodic function h{e), the Fourier coefficients of 
which occur explicitly in the expansion of V{e) stated in Thm. 5.1. 

Acknowledgements. The authors are grateful to Machiel van Frankenhuysen for 
his comments on a preliminary version of this paper. In particular, for indicating 
the current, more elegant and concise presentation of the main result (5.1). We 
also wish to thank Victor Shapiro for a helpful discussion concerning Fourier series, 
especially with regard to the discussion surrounding (4.8). 



Consider an approximation to the inner e-neighbourhood of the Koch curve, as 
shown in Fig. 3. Although we will eventually compute the neighborhood for the 
entire snowflake, we work with one third of it throughout the sequel (as depicted 
in the figure). 



T^dn — {D + inp : n e Z} U {inp : n G Z}. 



(1.11) 



2. Estimating the area 







Figure 2. The first four stages in the geometric construction of K . 
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Figure 3. An approximation to the inner e -neighbourhood of the Koch curve, with 
e £ l2- The refinement level here is based on the graph K2, the second stage in the 
geometric construction of the Koch curve (see Fig. 2). 



We will determine the area of the e-neighbourhood with functions that give the 
area of each kind of piece (rectangle, wedge, fringe, as shown in the figure) in 
terms of e, and functions that count the number of each of these pieces, in terms 
of e. As seen by comparing Fig. 3 to Fig. 4, the number of such pieces increases 
exponentially. 

We take the base of the Koch curve to have length 1 , and carry out our approx- 
imation for different values of e, as e — > 0"*". In particular, let 

/„ := (3-("+iV%/3,3-"/V3]. (2.1) 

Whenever e = 3^"/a/3, the approximation shifts to the next level of refinement. 
For example. Fig. 3 shows e £ I2, and Fig. 4 shows e £ I3. Consequently, for e G /o, 
it suffices to consider an e-neighbourhood of the prefractal curve Kq , and for e £ Ii, 
it suffices to consider an e-neighbourhood of the prefractal curve Ki, etc. Fig. 2 
shows these prefractal approximations. 

In general, we will discuss a neighbourhood of Kn, and we define the function 

1 



n{e) 



(2.2) 



to tell us for what n we have e £ In- Here, the square brackets indicate the floor 
function (integer part) and 

a; := -log3(e\/3), (2.3) 




Figure 4. A smaller e-neighbourhood of the Koch curve, fore £ I3. This refinement 
level is based on the graph K3, the third stage in the construction of the Koch curve. 
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a notation which will frequently prove convenient in the sequel. Further, we let 

{x} := X - [x] e [0, 1) (2.4) 

denote the fractional part of x. 

Observe that for e G /„ , n is fixed even as e is changing. To see how this is useful, 
consider that for all e in this interval, the number of rectangles (including those 
which overlap in the corners) is readily seen to be the fixed number 

r„ := 4". (2.5) 

Also, each of these rectangles has area e3"", where n is fixed as e traverses /„. 
Continuing in this constructive manner, we prove the following lemma. 

Lemma 2.1. For e € In, there are 
(i) "Tn = 4" rectangles, each with area £3~", 
(a) Wn = I (4" — 1) wedges, each with area 
(Hi) u„ = |(4" + 2) triangles, each with area and 
(iv) 4" components of fringe, each with area -^9"". 

Proof. We have already established (i). 

For (ii), we exploit the self-similarity of the Koch curve K to obtain the recurrence 
relation w„ = 4w„_i + 2, which we solve to find the number of wedges 

ri-l 

u;„:=^2.4^ = |(4"-l). (2.6) 

The area of each wedge is clearly 7r£^/6, as the angle is always fixed at 7r/3. 

(iii) To prevent double-counting, wc will need to keep track of the number of 
rectangles that overlap in the acute angles so that we may subtract the appropriate 
number of triangles 

n-l 

Un := 4" - ^ 4^ = I (4" + 2) . (2.7) 

Each of these triangles has area 

(iv) To measure the area of the fringe, note that the area under the entire Koch 
curve is given by •\/3/20, so the fringe of Kn will be this number scaled by (3~") . 
There are 4" components, one atop each rectangle (see Fig. 3). □ 

Lemma 2.1 gives a preliminary area formula for the £-neighbourhood. Here, 'pre- 
liminary' indicates the absence of the 'error estimate' developed in §3. 

Lemma 2.2. The e -neighbourhood of the Koch curve has approximate area 
V{e) = £2-^4-{-} (#9{-> + ^3{-> + i (l - Vs) ) - ^ (l + 2 Vs) . (2.8) 

This formula is approximate in the sense that it measures a region slightly larger 
than the actual £-neighbourhood. This discrepancy is accounted for and analyzed 
in detail in §3. 
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Proof of Lemma 2.2. Using (2.2) and (2.3), we obtain: 

4-=ie-^, 9-- = 3^2, (|)-=^^el-c, (|)- ^ 3^2-1?^ (2.9) 

Now using n = [x] = X — {x}, we compute the contributions of the rectangles, 
wedges, triangles, and fringe, respectively, as 

t/(£) = e(|)"^£2-D^^.4-W3M^ 



2 

V^{e) = ^(4" - 1) = . 4-W -s^^, 



14(e) = s!^ (4- + 2) = £2-^?^ . 4-{-} +s^^, and 

Putting all this together, V = Vr + Vyj — Vu + Vf gives the result. □ 

Remark 2.3. It is pleasing to find that this is in agreement with earlier pre- 
dictions of what V{e) should look like. In particular, [La-vFl, p. 209] gives the 
estimate 



V{e) « £2-13 A^4-M ggW ^ g . 3W _ 



(2.10) 



which differs only from our formula for V{e) in (2.8) by some constants and the 
final term of order e"^ . 

In (3.9) wc will require the Fourier scries of the periodic function e~^'^~^^V{e) 
as it is given by (2.8), so we recall the formula 

a-{^} = ^y ^ . (2.11) 

a ^-^ log a + 2'Kin 

This formula is valid for a > 0,a ^ 1 and has been used repeatedly in [La-vFl]. 
Note that it follows from Dirichlet's Theorem and thus holds in the sense of Fourier 
series. In particular, the series in (2.11) converges pointwise; this is also true in 
(2.13) and (2.14) below. 
We will make frequent use of the following identity in the sequel: 



g27rmx 



(£^3) '"^ = (-1)"£-"P, for n e Z, (2.12) 

where p = 2tt/ log 3 is the oscillatory period as in Thm. 1.1. With x = - log^ieVS), 
we can rewrite the Fourier expansion of a^^^^ given in (2.11) as 



a 



-M = ^V^^i^. (2.13) 
alog3^1og3a + znp 



Recall that {a;} = x — [x] denotes the fractional part of x. With D = logg 4, we 
use (2.13) to express (2.8) as a pointwise convergent Fourier series in e: 

^ ~ 31og3 \2^{D-2+inp) 2^{D-l+inp) ^ 2i'(D+mp) J ^ 
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3. Computing the error 

Now we must account for all the little 'trianglets', the small regions shaped like 
a crest of water on an ocean wave. These regions were included in our original 
calculation, but now must be subtracted. This error appeared in each of the rect- 
angles counted earlier, and so we refer to all the error from one of these rectangles 
as an 'error block'. Fig. 5 shows how this error is incurred and how it inherits a 
Cantoresque structure from the Koch curve. Actually, we will see in §3.2 that this 
'error' will have several terms, some of which are of the same order as the leading 
term in ^(e), which is proportional to e'^~^ by (2.14). Hence, caution should be 
exercised when carrying out such computations and one should not be too quick to 
set aside terms that appear negligible. 

3.1. Finding the area of an 'error block' 

In calculating the error, we begin by finding the area of one of these error blocks. 
Later, we will count how many of these error blocks there are, as a function of e. 
Note that n is fixed throughout §3.1, but s varies within /„. We define the function 

w = w(e) :=3-" = 3-[^l. (3.1) 

This function w{e) gives the width of one of the rectangles, as a function of e (see 
Fig. 6). Note that w{£) is constant as e traverses /„, as is n = n(e). From Fig. 6, 
one can work out that the area of both wedges adjacent to Ak is 



and that the area of the triangle above is 




Then for fc = 1, 2, . . . , the area of the trianglet Ak is given by 

. / N w(e) 9 . 1 / w(e) \ w(e) I ( wis) \ ^ ,„ „, 

and appears with multiphcity 2*^^^, as in Fig. 5. We use (3.1) and (2.3) to write 
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w{s) = 3-^3{^> = eVSi^y^^'h and define 

ox to o{x} — k+l/2 

"'fc SI'S ^ 

Hence the entire contribution of one error block may be written as 



fc=i 



,x N 2 ■ 



(3.3) 



(3.4) 



Recall the power series expansions 



sm 



-1 (2m)! 11^" + ^ , fr 

^=2^ 2^^(mV(2m+l) ^nd Vl 



-, _ \ ' (2m)!M'''"+^ 
^ 2^'"+im!(m+l)! ' 



which are valid for |u| < 1. We use these formulae with u = 2We' convergence 
is guaranteed by 

n < ^. — 3'"'V3 < VI ^ 1 

— 2We 2-3'' — 2 ^ ^' 

and the fact that the series in (3.4) starts with fc = 1. Then (3.4) becomes 



i?(£)=^2^-i 

oo 



(2m)! (3^)'^'"+^ 



1)! 24"'+i(m!p(2m+l) 
m— m— 



K _i_ \ " ^..fc; 
2 Z^ 2'i'"+-*m!(m+ 



fc=l 

oo oo 

EE 

oo 



(2m-2)!(3^)''"'+i _ (2m)! (3^)""'+^ 

/ . 2''™(m-l)!m! 2-4'" + i(m!)2(2m+l) 



(2m-2)! (3{=^} + i/2)2m+l / ^ (2m-l)2m \ ^2 
l^-*- 2m(2m+l) ^ 



(32m+i)fc 2'»'"(m-l)!m! 



E l / 1 A (2m-2)! (^3)""+^ f 1 2 
^(m — l)!m!(2m+l) V 3^"i + i / 

m— 1 

oo 

E (2™-2)! / I ^2 /o 

2^*^ ^ (m — l)!m !(2m+l)(3^'"+i-2) V3^" + i/ ^ " 



The interchange of sums is validated by checking absolute convergence of the final 
series via the ratio test, and then applying Fubini's Theorem to retrace our steps. 



>v(e)=3"" 

1 







A / 


1 




H 




V / 


\ 

\ 


1 













Figure 6. Finding the height of the central triangle. 
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3.2. Counting the error blocks 

Some blocks are present in their entirety as e traverses an interval /„ , while others 
are in the process of forming: two in each of the peaks and one at each end (see 
Fig. 7). Using the same notation as previously in Lemma 2.1, we count the complete 
and partial error blocks with 

c„ = r„ - u„ = i (4" - 4) , and p„ = u„ = § (4" + 2) . 

By means of (2.2)-(2.4), we convert c„ and Pn into functions of the continuous 
variable e, where e > 0: 

c(e) = i (^4-{"> -4) and p(e) = | (£^4-{"> + 2) . (3.6) 

With B{e) given by (3.5), the total error is thus 

Eie)^B{e)[c{e)+pie)hie)]. (3.7) 



Remark 3.1. The function h{e) in (3.7) is some periodic function that oscil- 
lates multiplicatively in a region bounded between and 1, indicating what portion 
of the partial error block has formed; see Fig. 5 and Fig. 7. We do not know h{e) 
explicitly, but we do know by the self-similarity of K that it has multiplicative 
period 3; i.e., h{e) = ft.(|). Using (2.12), the Fourier expansion 

He) = 5o(-l)"£-*"P = E = 9{x) (3.8) 

shows that we may also consider h{e) as an additively periodic function of the 
variable x = — \og^{eV5), with additive period 1. We refer the interested reader to 
§6 below for a further discussion of h{s), including a sketch of its graph, justification 
of the convergence of (3.8), and a brief discussion of some of its properties. 

We now return to the computation; substituting (3.5) and (3.6) into (3.7) gives 
E{e) = B{e) [^4-{"> (i + h{e)) + | {h{e) - 1) 




Figure 7. Error block formation. The ends are counted as partial because three of these 
pieces will be added together to make the entire snowflake. 
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C30 

1 \^ (2m-2)!(/t(e) + l/2) ( 4 2-D 

3 24'"-i(m-l)!m!(2m+l)(3^'"+i-2) ^ 

m=l 



, 1 (2m-2j!(/i(£)-l) / 1 ^2 

~^ 3 Z-> 24"--3(m-l)!m!(2m+l)(32-+i-2) Vs^^^^ ^ 
m=l 

00 

EV^ (2m-2)l(4-:i^'" + i)(-l)"(ft(£) + l/2) ^2-g-mp 

/ . 2*™+i(m-l)!m!(2ni+l)(3^'" + '--2)(L>-2in-l+mp)^ 



31 _ 

ni=l rteZ 

CX) 

^_ V- V- (2m-2)!(l-3^"'+^)(-l)"(fc(^)-l) ^2-mp 

^ 3 log 3 2*"'-='(m-l)!m!(2m+l)(3^'"+i-2)(-2m-l+mp)'^ 

m=l nGZ 



+ TO E('^(^) - 1)(-Tn)(-1)"£'-'"^ (3.9) 

where we have written the constants 6„ and r„ in the shorthand notation as foUows: 
00 

. \^ (2m-2)!(3^'"+^-4) 

■~ 2*"'+i(TO-l)!m!(2m+l)(3^'"+i-2)(D-2TO-l+mp)' l,o.lu; 

m=l 

00 

, _ (2m-2)!(3'''"+^-i) 

ana r„ :— 24™-3(TO-l)!m!(2m+l)(32'"+i-2)(-2m-l+mp) • K'^-'-'-) 



In the third equahty of (3.9), we have apphed (2.13) to a = 4/3^™+^ and to 
a = 1/3^™+^, respectively. By the ratio test, the complex numbers 6„ and t„ given 
by (3.10) and (3.11) arc well-defined. This fact, combined with Fubini's Theorem, 
enables us to justify the interchange of sums in the last equality of (3.9). 



4. Computing the area 

Now that we have estimate (2.14) for the area of the neighbourhood of the Koch 
curve V{s), and formula (3.9) for the 'error' E{e), we can find the exact area of the 
inner neighbourhood of the full Koch snowflake as follows: 

V{6) = 3 (V{6) - E{s)) 

= ^ V ( - -^'^^ -4- , 7r-3'V^ \ -,n-2-D-mp 

log 3 Z^ (25(r>-2+mp) 2'-'{D-l+inp) ^ W{D+inp) J ^ ' ^ 

+ I3i3 E(^(^) + l/2)(-l)"6„£^-^-»P 

J + 2V3) e^ + j^, Y.^h{e) - 1)(-I)"r„e2--P 

raeZ 

\2^(D-2+in^) 23(£)-l+mp) 2a(D+mp) 

rteZ 

+ % + /i(£)6„) (_i)V-"P 
+ Ei3 E (-^" - log 3 (5 + '^0 + '^(^)^") (-l)"^'""^ (4-1) 



nez 
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where 6q is the Kronecker delta. Therefore, we have 

y(£) = Gl(e)£2-^ + G2(£)£^ (4.2) 
where the periodic functions Gi and G2 are given by 

Gi{e) : = ^ (a„ + 6„ft(£)) (-1)"£-»p (4.3) 

and G2(e) : = E + ^"'^(^)) (-l)"^"'"''- (4-4) 

Here wc have used (3.10) and (3.11), and introduced the notation a„ and (t„ (stated 
exphcitly in (5.3)). We wish to rearrange these series so as to collect all factors of 
s. First, we can split the sum in (4.3) as 

^l(^) = E «n(-l)"£-^"P + ^ E ^n(-l)"£-^"P (4.5) 

because a{s) := '^"("l)"^"*"^ ^^'^ K^) •= Snez ^"("1)"^"*"'' 

convergent: a(e) converges for the same reason as (2.14), and one can show that 

00 

- V (41og3)(2m-2)! ( 4 . ^'W (A (i\ 

— 2^ 24™+i3™-i/2(m_i)!m!(2m+l)(32™+i-2) \W^^ ) V^'"^ 

converges to a well-defined distribution induced by a locally integrable function; 
one proves directly that < c/|n| by writing (3.10) as 

oo oo 
^" = E D-2nfri+,„p , With E /3„ < OO. (4.7) 

m=l m— 1 

Then the rearrangement leading to (4.6) is justified via the "descent method" with 
g = 2, as described in Rem. 4.1 below. Note that the right-hand side of (4.6) 
converges by the ratio test and is thus defined pointwise on M ~ Z. Therefore one 
also sees that 6(e) and h{e) are periodic functions. Considered as functions of the 
variable x = log3(l/eA/3), both have period 1 with h continuous for < a; < 1 and 
h continuous for < .t < 1. Further, each is monotonic on its period interval, and 
possesses a bounded jump discontinuity only at the endpoint. This may be seen for 
6(e) from (4.6) and for hie) from §6. 

Although we have used distributional arguments to obtain (4.6), the right-hand 
side of (4.6) is locally integrable and has a representation as a piecewise continuous 
function. Thus the distributional equality in (4.G) actually holds pointwise in the 
sense of Fourier series and all our results are still valid pointwise. Recall that the 
Dirichlet- Jordan Theorem [Zy, Thm. II. 8.1] states that if / is periodic and (locally) 
of bounded variation, its Fourier series converges pointwise to (/(.t— ) + /(.t+)) /2. 
As described above, 6(e) and hie) are each of bounded variation and therefore 
by [Zy, Thm. II.4.12] we have 6„ = 0(l/n) and Qn = 0{l/n) as n — > ±oo. Finally, 
[Zy, Thm. IX.4.11] may be applied to yield the pointwise equality 

b{e)h{e) = J2Y1 &a5n-a(-l)"£-'"P. (4-8) 

neZctez 

This theorem applies because (4.6) shows that 6(e) is bounded away from 0. 
Now that all the e's are combined, we substitute (4.8) back into (4.5) and rewrite 
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Gi{e) = {a{e) + b{e)h{e)) /log 3 as in (5.2a): 




— inp 



(4.9) 



Manipulating G2{s) similarly, we are able to rewrite (4.4) in its final form (5.2b), 
and thereby complete the proof of Thm. 5.1. 

Remark 4.1. The convergence of the Fourier series associated to a periodic 
distribution is proved via the descent method by integrating both sides q times 
(for sufficiently large q) so that one has pointwise convergence. After enough inte- 
grations, the distribution will be a smooth function, the series involved will con- 
verge absolutely, and the Weierstrass theorem can be applied pointwise. At this 
point, rearrangements or interchanges of series are justified pointwise, and wc ob- 
tain a pointwise formula for the q**^ antiderivative of the desired function. Then one 
takes the distributional derivative q times to obtain the desired formula. Sec [La- 
vFl, Rem. 4.14]. How large the positive integer q needs to be depends on the order 
of polynomial growth of the Fourier coefficients. Recall that the Fourier series of a 
periodic distribution converges distributionally if and only if the Fourier coefficients 
arc of slow growth, i.e., do not grow faster than polynomially. Moreover, from the 
point of view of distributions, there is no distinction to be made between convergent 
trigonometric series and Fourier series. See [Sch, §VII,I]. 



We can now state our main result in the following more precise form of Thm. 1.1: 

Theorem 5.1. The area of the inner e -neighbourhood of the Koch snowflake is 
given pointwise by the following tube formula: 



where Gi and G2 are periodic functions of multiplicative period 3, given by 



5. Main results 



y(£) = Gl(£)e2-^ + G2(e)£^ 



(5.1) 




(5.2b) 



(5.2a) 



where an,b„,a„, and t„ are the complex numbers given by 

35/2 33/2 7r~33/2 1 



2^{D-2 + inp) 23{D-l + inp) 23{D + inp) 2 




(5.3) 
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Figure 8. The possible complex dimensions of K and dO.. The Minkowski dimension is 
D = logg 4 and the oscillatory period is p — . 



where Sq = 1 and Sq = for n ^ is the Kronecker delta. 

In Thm. 5.1, D = log3 4 is the Minkowski dimension of the Koch snowflake dVl and 
p = 27r/log3 is its oscillatory period, following the terminology of [La-vFl]. The 
numbers ga appearing in (5.2) are the Fourier coefficients of the periodic function 
h{e), a suitable nonlinear analogue of the Cantor-Lebesgue function, defined in 
Rem. 3.1 and further discussed in §6 below. 

The reader may easily check that formula (5.1) can also be written 

V{e) = VnS^^''''^'' + E ^"^'^"^ (5-4) 

for suitable constants (pmi^n which depend only on n. Now by analogy with the 
tube formula (1.8) from [La-vFl], we interpret the exponents of e in (5.4) as the 
'complex co-dimensions' of dfl. Hence, we can simply read off the possible complex 
dimensions, and as depicted in Fig. 8, we obtain a set of possible complex dimensions 

'Don = {D + inp : n e Z} U {inp : neZ}. (5.5) 

One caveat should be mentioned: this is assuming that none of the coefficients ipn 
or i/jn vanishes in (5.4). Indeed, in that case we would only be able to say that the 
complex dimensions are a subset of the right-hand side of (5.5). We expect that, 
following the 'approximate tube formula' obtained in [La-vFl], the set of complex 
dimensions should contain all numbers of the form D + inp. 

Remark 5.2. We expect our methods to work for all lattice self-similar frac- 
tals (i.e., those for which the underlying scaling ratios are rationally dependent) as 
well as for other examples considered in [La-vFl] including the Cantor-Lebesgue 
curve, a self-afhne fractal. For example, we have already obtained the counterpart of 
Thm. 5.1 for the square (rather than triangular) snowflake curve. By applying den- 
sity arguments (as in [La-vFl, Chap. 2]), our methods may also yield information 
about the complex dimensions of nonlattice fractals. 



A TUBE FORMULA FOR THE KOCH SNOWFLAKE CURVE. 



15 



Remark 5.3. It should be pointed out that in the present paper, we do not 
provide a direct definition of the complex dimensions of the Koch snowflake curve 
dCl (or of other fractals in M^). Instead, we reason by analogy with formula (1.8) 
above to deduce from our tube formula (5.1) the possible complex dimensions of dfl 
(or oi K). As is seen in the proof of Thm. 5.1, the tube formula for the Koch curve K 
is of the same form as for that of the snowflake curve dil. It follows that K and dil 
have the same possible complex dimensions (see Cor. 1.2). In subsequent work, we 
plan to define the geometric zeta function (qq = Cdn{s) in this context to actually 
deduce the complex dimensions of dil directly as the poles of the meromorphic 
continuation of Can- 

In fact, this project is well underway. A self-similar tiling defined in terms of 
an iterated function system on M'' is constructed in [Pe]. In [LaPe2], we use this 
tiling and aspects of geometric measure theory discussed in [LaPel] to define a zeta 
function whose poles give the complex dimensions directly. We then use the zeta 
function and complex dimensions to obtain an explicit distributional inner tube 
formula for self-similar systems in M"^, analogous to [La-vFl, Chap. 6]. We stress 
that although Thm. 5.1 was helpful in developing [Pe] and [LaPel-2], it is not a 
consequence or corollary of this more recent work. 

Remark 5.4. In the long-term, by analogy with Hermann Wcyl's tube formula 
[We] for smooth Riemannian submanifolds (see [Gr]), we would like to interpret 
the coefficients a„,6„,(j„,r„ of the tube formula (5.1) in terms of an appropriate 
substitute of the 'Weyl curvatures' in this context. See the corresponding discussion 
in [La-vFl, §6.1.1 and §10.5] for fractal strings; also see [La-vF3]. 

This is a very difficult open problem and is still far from being resolved, even in 
the one-dimensional case of fractal strings. See [La-vFl, §6.1.1] and [LaPel-2]. 

Remark 5.5. (Reality principle.) As is the case for the complex dimensions 
of self-similar strings (see [La-vFl, Chap. 2] and [La-vF2-3]), the possible complex 
dimensions ofdCl come in complex conjugate pairs, with attached complex conjugate 
coefficients. Indeed, since = g-a (see §6), a simple inspection of the formulas in 
(5.3) shows that for every n € N, 

a:^ = a-n,bn = b-njOn = CT-n, a,ndT^ = T-n. (5.6) 

It follows that ao, 5o, (To, and tq are reals and that Gi and G2 in (5.2) are reai- valued, 
in agreement with the fact that V{e) represents an area. 

6. The Cantor-like function h{e). 
We close this paper by further discussing the nonlinear Cantor-like function 

h{e) = J2 ffa(-l)"^-^"P = E 5"^'""' = ai^)' (6.1) 

introduced in (3.8). Note that since h is real- valued (in fact, < ft. < /x < 1), we 
have g-a = g^ for all a £ Z. Further, recall from §3.2 that in view of the self- 
similarity of K, h{e) is multiplicatively periodic with period 3, i.e., h{e) = /i(|). 
Alternatively, (6.1) shows that it can be thought of as an additively periodic function 
of X with period 1, i.e., g{x) = g{x + 1). By the geometric definition of h{e), we 
see that it is continuous and even monotonic when restricted to one of its period 
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C(e) 

B{e) _ 

Figure 9. /i is the ratio vol2(C(e))/ vol2(-B(e)). 

intervals /„ := ^3-n-3/2 2-»i-i/2j ^ gi^^e h is of bounded variation, its Fourier 
series converges pointwise by [Zy, Thm. II. 8.1] and its Fourier coefficients satisfy 

5a = 0(1/q;), as a ^ ±oo (6.2) 

by [Zy, Thm. II. 4. 12], as discussed in §4. 

Further, since h{£) is defined as a ratio of areas (see Remark 3.1), we have h{£) G 
[0, ijl) for all e > 0. Note that /i(e) < /i < 1 and so /i(e) does not attain the value 
1; the error blocks being formed are never complete. The partial error blocks only 
form across the first | of the line segment beneath them. They reach this point 
precisely when e = 3~"/-y3 for some n > 1. Back in (3.4), we found the error of a 
single error block to be given by 



i?(£)=^2^-lA,(£), 



fc=l 



where Ak{e) is given by (3.2). Thus the supremum of h{e) will be the ratio 

which will be the same constant for each eu = 3^'^^^''^, fc = 1, 2, . . . (see Fig. 9). 
In other words, the number we need is 



Ai{et) + \T.Z2'^'-'Meu) Ai(£,) + Er^,2^-^Afc(£fc) 



Note that although this definition of initially appears to depend on fc, the ratio in 
question is between two areas which have exactly the same proportion at each e^; 
this is a direct consequence of the self-similarity of the Koch curve. In other words, 
if C(e) is the area indicated in Fig. 9, then the relations 

3C(f ) = C{e) and 3B(|) = B{e) 

show that fi is well-defined. 

We now approximate hie) by a function which shares its essential properties: 

(i) h{ek) = lim h{ek + 'd)^0, 
(u) hm h(ek + ■&) = u, 

where for fc = 1, 2, . . . , = again. That is, h{e) goes from to as e goes 
from 3"*^ to 3"^'^+^^. Using again the notation x ~ — log;^{ey/3) and {a;} — x — [x], 
we see that the function 

hie) ^ ^1 ■ {-[x] - x} (6.3) 

shares both of these properties but is much smoother. Indeed, h{e) only has points 
of nondifferentiability at each and is otherwise a smooth logarithmic curve; see 
Fig. 10. The true h(e), by contrast, is a much more complex object that deserves 
further study in later work. 
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